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In  our  recent  work,  a  multi-physics  computational  model  for  the  conventional  gas  metal  arc  welding 
(GMAW)  joining  process  was  introduced.  The  model  is  of  a  modular  type  and  comprises  five  modules,  each 
designed  to  handle  a  specific  aspect  of  the  GMAW  process,  i.e.:  (i)  electro-dynamics  of  the  welding-gun;  (ii) 
radiation-/convection-controlled  heat  transfer  from  the  electric-arc  to  the  workpiece  and  mass  transfer 
from  the  filler-metal  consumable  electrode  to  the  weld;  (iii)  prediction  of  the  temporal  evolution  and  the 
spatial  distribution  of  thermal  and  mechanical  fields  within  the  weld  region  during  the  GMAW  joining 
process;  (iv)  the  resulting  temporal  evolution  and  spatial  distribution  of  the  material  microstructure 
throughout  the  weld  region;  and  (v)  spatial  distribution  of  the  as-welded  material  mechanical  properties.  In 
the  present  work,  the  GMAW  process  model  has  been  upgraded  with  respect  to  its  predictive  capabilities 
regarding  the  spatial  distribution  of  the  mechanical  properties  controlling  the  ballistic-limit  (i.e.,  pene¬ 
tration-resistance)  of  the  weld.  The  model  is  upgraded  through  the  introduction  of  the  sixth  module  in  the 
present  work  in  recognition  of  the  fact  that  in  thick  steel  GMAW  weldments,  the  overall  ballistic  perfor¬ 
mance  of  the  armor  may  become  controlled  by  the  (often  inferior)  ballistic  limits  of  its  weld  (fusion  and 
heat-affected)  zones.  To  demonstrate  the  utility  of  the  upgraded  GMAW  process  model,  it  is  next  applied  to 
the  case  of  butt-welding  of  a  prototypical  high-hardness  armor-grade  martensitic  steel,  MIL  A46100.  The 
model  predictions  concerning  the  spatial  distribution  of  the  material  microstructure  and  ballistic-limit¬ 
controlling  mechanical  properties  within  the  MIL  A46100  butt-weld  are  found  to  be  consistent  with  prior 
observations  and  general  expectations. 


Keywords  ballistic  limit,  gas  metal  arc  welding  (GMAW)  process 
modeling,  welded  all-metal  armor 


1.  Introduction 

Our  recently  proposed  multi-physics  computational  model  for 
the  conventional  gas  metal  arc  welding  (GMAW)  joining  process 
(Ref  1  -3),  Fig.  1 ,  is  extended  in  the  present  work  in  order  to  enable 
predictions  of  the  spatial  distribution,  throughout  the  weld,  of  the 
mechanical  properties  controlling  the  ballistic-limit.  The  welding 
process  parameters  were  already  correlated  with  the  spatial 
distribution  of  the  material  microstructure  and  properties  within 
the  weld  (consisting  of  the  solidified  weld  pool,  also  referred  to  as 
the  fusion  zone,  FZ,  and  the  adjacent  heat-affected  zone,  HAZ) 
using  the  original  GMAW  process  model.  However,  regarding  the 
prediction  of  the  spatial  distribution  of  the  mechanical  properties 
that  control  the  ballistic  limit  in  thick  metallic-armor  weldments, 
the  model  had  serious  limitations.  Specifically,  the  model  was 
unable  to  predict  spatial  distribution,  throughout  the  weld,  of  the 
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following  material  properties:  (a)  material  strength,  (b)  stiffness, 
(c)  volume-based  specific  heat,  (d)  exponent  of  strain-hardening, 
(e)  exponent  of  strain-rate  sensitivity,  and  (f)  material  shear- 
strength  temperature-dependence.  In  our  previous  work  (Ref  1-3), 
the  model  was: 

(a)  applied  to  the  GMAW  of  MIL  A46100  (an  air-quenchable, 
self-tempered,  high-hardness,  low-alloy  martensitic,  ar¬ 
mor-grade  steel).  A  summary  of  Johnson-Cook  material- 
model  parameters  and  other  thermo-mechanical  properties 
for  this  material  is  given  in  Tables  1  and  2;  and 

(b)  validated  by  comparing  the  model  predictions  (regarding 
the  spatial  distribution  of  the  material  micro  structure 
throughout  the  weld)  with  their  experimental  counterparts. 
In  order  to  predict  spatial  distribution  of  the  mechanical 
properties  controlling  the  ballistic  limit  in  welded-armor 
of  this  material,  the  extended  GMAW  process  model  is 
applied  to  the  same  steel  grade  in  the  present  work. 

The  main  objective  of  the  present  work  is  to  improve  the 
capability  of  our  recently  proposed  GMAW  process  model  to  predict 
the  spatial  distribution  of  material  properties  within  the  weld  region. 
This  capability  is  important  for  two  reasons.  First,  these  properties 
may  often  be  inferior  relative  to  their  base-metal  counterparts. 
Second,  the  width  of  the  weld  in  thick  metallic-armor  is  often 
comparable  to  the  armor  thickness,  and  therefore  may  represent  a 
significant  portion  of  the  armor  exposed-surface  area  (Ref  4). 

A  concise  summary  of  the  multi-physics  GMAW  process 
model  and  its  five  modules,  as  reported  in  Ref  1-3,  is  presented 
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Fig.  1  A  schematic  representation  of  the  conventional  gas  metal  arc  welding  (GMAW)  process 


Table  1  Johnson-Cook  strength  model  material  parameters  for  MIL  A46100 


Parameter 

Symbol 

Units 

Value 

Young’s  modulus 

E 

GPa 

205-215 

Poisson’s  ratio 

V 

N/A 

0.285-0.295 

Reference  strength 

A 

MPa 

1000-1100 

Strain-hardening  parameter 

B 

MPa 

250.0 

Strain-hardening  exponent 

n 

N/A 

0.12 

Strain-rate  coefficient 

C 

N/A 

0.02 

Room  temperature 

T 

1  room 

K 

298.0 

Melting  temperature 

Tnelt 

K 

1720 

Temperature  exponent 

m 

N/A 

0.5 

Table  2  General,  thermo-mechanical,  and  thermal  parameters  for  MIL  A46100 

Parameter 

Symbol 

Units 

Value 

Material  mass  density 

P 

kg/m3 

7840-7860 

Coefficient  of  linear  thermal  expansion 

a 

1/K 

lle-6-12e-6 

Specific  heat 

Cp 

J/kg  K 

440-520 

Thermal  conductivity 

k 

W/mK 

35-50 

Heat  transfer  coefficient 

h 

W/m2  K 

45 

Sink  temperature 

Link 

K 

298 

Emissivity 

£ 

N/A 

0.57 

Ambient  temperature 

T 

1  amb 

K 

298 

in  section  2.  In  section  3.1,  a  brief  overview  is  provided  of  the 
major  changes  made  in  the  fifth  module,  the  microstructure/ 
property  relationship  module.  Details  related  to  construction  of 
the  sixth  module,  the  ballistic-limit  prediction  module  is 
presented  in  section  3.2.  In  the  same  section,  the  results 
pertaining  to  the  spatial  distribution  of  the  ballistic  limit 
throughout  the  weld  in  MIL  A46100  steel  are  next  presented 
and  discussed.  The  main  conclusions  resulting  from  the  present 
work  are  summarized  in  section  4. 

2.  Multi-physics  GMAW  Process  Model 

In  order  to  help  clarify  the  structure  of  the  present  six- 
module  GMAW  process  model,  a  flowchart  revealing  the 


sequence  and  interconnectivity  of  the  modules  is  depicted  in 
Fig.  2.  The  first  five  modules  were  developed  in  Ref  1-3  and, 
hence,  no  summary  of  the  key  physical  concepts  and  principles 
pertaining  to  these  modules  will  be  provided  here.  Instead,  a 
few  prototypical  results  yielded  by  these  modules  will  be 
presented  and  discussed.  The  sixth  module  shown  in  Fig.  2  has 
been  developed  in  the  present  work. 

2.1  Welding-Gun  Module 

Figure  3(a)  and  (b)  shows  typical  results  of  the  temporal 
evolutions  of  the  welding  voltage  (input)  and  the  welding 
power  (output),  as  yielded  by  this  module.  It  should  be  noted 
that  the  output  welding  power  shows  considerable  fluctuations 
even  when  the  input  welding  voltage  is  set  to  a  constant  value. 
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Input 


Welding-Gun  Module  (1) 


•  Open-circuit  Voltage;  Electrode-tip  /  Weld  Distance;  Filler-metal  Feed  Speed; 
Gun  Travel  Speed 

Output 

•  Heat  Power;  Molten  Filler-metal  Deposition  Rate 


Fig.  2  The  structure  of  the  multi-physics  GMAW  process  model  consisting  of  six  modules 


Consequently,  this  module  establishes  the  functional  relation¬ 
ship  between  the  input  voltage  and  mean  output  power,  and  this 
functional  relationship  is  subsequently  passed  to  the  next 
module. 

2.2  Electric- Arc- Based  Heat-Source/Metal-Spray  Module 

Figure  4  depicts  a  three-dimensional  surface  plot  of  nor¬ 
malized  surface  heat-flux  as  a  function  of  distances  along  the 
weld  transverse  direction  (x)  and  the  weld  longitudinal  direction 
(y).  It  should  be  noted  that  x  and  y  (as  well  as  the  standard- 
deviation  of  the  assumed  normal-circular  power-density  distri¬ 
bution  function)  are  all  normalized  by  the  (equal  transverse  and 
longitudinal)  workpiece  half-edge  lengths. 


2.3  Thermo-mechanical  GMAW  Process  Module 

The  results  presented  below  were  obtained  for  the  case  of 
MIL  A46100  GMAW  butt-welding. 

2.3.1  Temporal  Evolution  of  the  Weldment  Temperature 
Field.  Typical  results  pertaining  to  the  temporal  evolution  of 
the  temperature  field  within  the  weld  region  over  the  (explicitly 
analyzed)  right  half  of  the  weldment  are  depicted  in  Fig.  5(a)- 
(d),  along  with  the  coordinate  system  used  in  this  analysis.  The 
results  shown  in  Fig.  5(a)-(d)  are  obtained  at  relative  welding 
times  of  0.6,  2.1,  3.6,  and  4.8  s,  respectively.  The  GMAW 
process  parameters  used  to  obtain  these  results  are  as  follows: 
welding  input  voltage  =  30  V,  welding  current  =  200  A,  elec¬ 
trode  diameter  =  1  mm,  electrode-tip/weld  distance  =1.3  cm, 
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Fig.  3  Typical  results  yielded  by  the  welding-gun  module  pertain¬ 
ing  to  temporal  evolutions  of  (a)  the  input  voltage;  and  (b)  the  out¬ 
put  power 


electrode  feed-rate  =10  cm/s,  and  gun  travel  speed  =  1  cm/s. 
It  should  be  noted  that  regions  of  the  weldment  with  a 
temperature  exceeding  the  liquidus  temperature  (defined  as  the 
temperature  at  which  the  volume  fraction  of  the  liquid  phase, 
during  cooling  from  the  melt,  first  begins  to  deviate  from 
100%)  are  displayed  as  red,  for  improved  clarity. 

Examination  of  the  results  displayed  in  Fig.  5(a)-(d)  reveals 
that:  (a)  after  a  brief  transient  period,  the  FZ  acquires  a  nearly 
constant  size  and  shape  as  it  travels  along  the  welding  direction 
(following  the  position  of  the  weld  gun);  (b)  as  welding 
proceeds,  the  previously  molten  material  within  the  FZ 
solidifies  (and  continues  to  cool)  due  to  natural  convection 
and  radiation  to  the  surroundings,  together  with  conduction 
through  the  adjacent  workpiece  material  region;  and  (c)  under 
the  aforementioned  welding  conditions,  the  FZ  extends  down¬ 
ward  into  the  workpiece  thickness  by  ca.  40-45%. 


2.3.2  Temporal  Evolution  of  Temperature  Within  the 
FZ  and  HAZ.  Figure  6(a)  and  (b)  shows  typical  results 
pertaining  to  the  temporal  evolution  of  temperature  at  fixed 
locations  within  the  FZ  and  HAZ,  respectively.  Results 
pertaining  to  the  weldment  (through- the-thickness)  mid-plane 
within  the  FZ  are  displayed  in  Fig.  6(a).  In  this  figure,  the 
distance  of  the  material  point  in  question  from  the  weld  y-z 
symmetry  plane  is  denoted  via  the  curve  labels.  Results 
pertaining  to  the  weldment  mid-plane  within  the  HAZ,  on  the 
other  hand,  are  displayed  in  Fig.  6(b).  The  curve  labels  used  in 
this  figure  represent  the  distance  of  the  material  point  in 
question  from  the  FZ/HAZ  interface.  It  should  be  noted  that 
before  the  temperature-history  results  generated  in  this  module 
are  passed  to  the  fourth  (microstructure-evolution)  module, 
they  are  computed  for  a  larger  number  of  closely  spaced 
material  points  in  the  FZ  and  HAZ. 

Examination  of  the  results  displayed  in  Fig.  6(a)  and  (b) 
reveals  that: 

(a)  Temperatures  in  excess  of  the  MIF  A46100  liquidus 
temperature  (1772  K),  are  experienced  by  material 
points  which  reside  within  the  FZ.  As  a  result  of  the  ar¬ 
rival  of  the  weld  gun,  these  points  are  subjected  to  a 
steep  rise  in  temperature,  which  is  followed  by  a  rela¬ 
tively  steep  drop  in  temperature  due  to  the  subsequent 
departure  of  the  weld  gun;  and 

(b)  Temperatures  which  are  in  excess  of  the  MIF  A46100 
Ac\  (defined  as  the  highest  temperature  at  which  austen¬ 
ite  is  still  present  during  cooling)  temperature  (982  K), 
but  which  are  lower  than  the  MIF  A46100  liquidus  tem¬ 
perature,  are  experienced  by  material  points  which  re¬ 
side  within  the  HAZ.  These  points  spend  more  time  at 
temperatures  near  the  peak  temperature,  are  subjected  to 
a  significantly  lower  heating  rate  and  undergo  relatively 
slow  cooling. 

2.3.3  Thermal-Strain/Residual-Stress  Fields.  High  ther¬ 
mal  stresses  and,  in  turn,  plastic  deformation  (permanent 
distortion)  can  occur  in  the  weld  region  if  non-uniformity  in  the 
temperature  distribution  throughout  this  region  is  sufficiently 
large.  In  particular,  due  to  highly  non-uniform  input  thermal 
flux  (as  seen  in  Fig.  4),  the  FZ  may  experience  high  thermal 
stresses  that  lead  to  eventual  plastic  deformation.  This  would 
cause  weldments  to  acquire  residual  stresses  upon  cooling  to 
room  temperature.  As  explained  earlier,  the  present  multi¬ 
physics  GMAW  process  model,  due  to  its  thermo-mechanical 
character,  is  capable  of  predicting  the  development  of  such 
welding-induced  permanent  distortions  and  residual  stresses. 
Figure  7(a)  and  (b)  shows  typical  thermal-strain/residual-stress 
results  obtained  using  the  present  rendition  of  the  thermo¬ 
mechanical  GMAW  process  module  (Ref  5-13).  Spatial  distri¬ 
bution  of  the  equivalent  plastic  strain  and  the  von  Mises 
equivalent  residual  stress,  over  a  transverse  section  of  the  weld, 
are  depicted  in  Fig.  7(a)  and  (b),  respectively.  Instead  of 
displaying  the  contour-level  legends,  the  maximum  values  of 
the  equivalent  plastic  strain  and  the  von  Mises  equivalent 
residual  stress  are  denoted  in  these  figures,  for  improved  clarity. 
Examination  of  the  results  displayed  in  Fig.  7(a)-(b)  show  that 
the  largest  plasticity/residual-stress  effects  are  observed  within 
the  FZ  and  the  adjacent  portions  of  the  HAZ,  as  is  to  be 
expected. 

Welding-induced  permanent  distortions  may  affect  func¬ 
tional  performance  of  a  weldment,  due  to  the  associated  shape 
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Fig.  4  A  normalized  heat-flux  three-dimensional  surface  plot  (for  the  case  of  the  circular-normal  distribution)  as  a  function  of  the  lateral  and 
longitudinal  locations  (normalized  by  the  corresponding  workpiece  half-edge  lengths)  on  the  workpiece  top  surface.  The  x  =  y  =  0  location,  asso¬ 
ciated  with  the  maximum  heat-flux,  corresponds  to  the  current  position  of  the  electrode  axis 


and  dimensional  changes,  while  the  presence  of  the  (tensile) 
residual  stresses  may  degrade  structural  performance  and/or 
affect  reliability/durability  of  the  weld.  In  addition,  through  the 
so-called  “deformation-induced”  and  “stress-assisted”  trans¬ 
formation  effects  (e.g.,  Ref  14),  plastic  strains  and  residual 
stresses  may  affect  the  progress  of  the  austenite  — »  bainite  and 
austenite  —>  martensite  displacive  phase  transformations. 
These  phase  transformations  will  be  presented  in  greater  detail 
in  the  section  dealing  with  the  microstructure-evolution  (fourth) 
module.  These  effects  are  not  included  in  the  current  rendition 
of  the  microstructure-evolution  module,  even  though  they  may 
somewhat  influence  the  predictions  of  the  micro  structure 
distribution  within  the  weld.  However,  as  demonstrated  in 
Ref  1 ,  within  the  common  GMAW  process  parameter  envelope, 
the  role  of  these  effects  is  relatively  small  and,  hence,  this  is 
presently  not  considered  as  a  major  issue  in  the  current  module. 

2.4  Microstructure-Evolution  Module 

Figure  8(a)-(f)  depicts  typical  results  yielded  by  the  micro- 
structure-evolution  module,  for  MIL  A46100  that  is  initially  in 
the  as-hot-rolled  and  self-tempered  martensitic  state.  These 
results  are  obtained  under  the  following  GMAW  process 
conditions:  input  welding  voltage  =  30  V,  welding  cur¬ 
rent  =  200  A,  electrode  diameter  =  1  mm,  electrode-tip/weld 
distance  =1.3  cm,  electrode  feed-rate  =10  cm/s,  and  gun 
travel  speed  =  1  cm/s.  Figure  8(a)-(e)  shows  the  spatial  distri¬ 
butions,  over  a  transverse  section  of  the  weld  and  the  workpiece 
region  adjacent  to  the  weld,  of  the  volume  fractions  for  the 
following  phases:  (a)  allotriomorphic  ferrite,  (b)  Widmanstatten 
ferrite,  (c)  bainite,  (d)  freshly  formed  martensite,  and  (e) 


tempered  martensite,  respectively.  The  corresponding  spatial 
distribution  of  the  prior-austenite  grain  size  is  shown  in 
Fig.  8(f). 

Examination  of  the  results  displayed  in  Fig.  8(a)-(f)  reveals 
that: 

(a)  The  FZ  consists  only  of  allotriomorphic  ferrite  and  Wid¬ 
manstatten  ferrite,  Fig.  8(a)-(b).  This  observation  is  to 
be  expected,  considering  the  fact  that  the  material  within 
this  region  is  initially  subjected  to  temperatures  signifi¬ 
cantly  higher  than  the  material  liquidus  temperature  and, 
hence,  is  associated  with  a  long  cooling  time.  Further¬ 
more,  since  austenite  chemical  composition  is  nearly 
identical  to  that  of  the  alloy  itself,  the  material  within 
the  FZ  possesses  relatively  low  hardenability  (the  ease 
of  transformation  of  austenite  into  martensite); 

(b)  There  is  a  fairly  small  (less  than  30  vol.%)  volume  frac¬ 
tion  of  freshly  formed  martensite  within  the  FZ, 
Fig.  8(d).  This  finding  is  fully  consistent  with  the  fact 
that  the  weld  region  experiences  relatively  low  cooling 
rates  and  contains  austenite  with  relatively  low  hardena¬ 
bility.  Upon  crossing  the  FZ/HAZ  interface  and  entering 
the  HAZ,  the  volume  fraction  of  freshly  formed  mar¬ 
tensite  first  increases  (to  a  value  in  excess  of  85  vol.%) 
and  then  decreases  (to  a  value  of  ca.  10  vol.%).  The 
volume  fraction  of  the  freshly  formed  martensite  subse¬ 
quently  drops  to  zero  upon  crossing  the  HAZ/unaf- 
fected-material  boundary.  A  compromise  between  the 
amount  of  austenite  available  to  transform  into  martens¬ 
ite  and  the  hardenability  of  austenite  is  responsible  for 
the  spatial  distribution  of  the  volume  fraction  of  freshly 
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Fig.  5  Spatial  distribution  of  the  temperature  field  in  the  weld  region  over  the  (analyzed)  right-portion  of  the  MIL  A46100  weldment  at  weld¬ 
ing  times  of  (a)  0.6  s,  (b)  2.1  s,  (c)  3.6  s,  and  (d)  4.8  s  under  the  following  welding  conditions:  welding  input  voltage  =  30  V,  welding  cur¬ 
rent  =  200  A,  electrode  diameter  =  1  mm,  electrode-tip/weld  distance  =  1.3  cm,  electrode  feed-rate  =10  cm/s,  and  gun  travel  speed  =  1  cm/s 
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Fig.  6  Temporal  evolution  of  temperature:  (a)  along  the  workpiece 
mid-plane  within  the  FZ  (the  curve  labels  denote  the  distance  of  the 
subject  material  point  from  the  weld  y-z  symmetry  plane);  and  (b) 
along  the  workpiece  mid-plane  within  the  HAZ  (the  curve  labels  de¬ 
note  the  distance  of  the  subject  material  point  from  the  HAZ/FZ 
interface) 

formed  martensite.  Specifically,  it  is  assumed  that  full 
austenitization  takes  place  in  the  HAZ  regions  which 
were  exposed  to  temperatures  above  Ac3  (defined  as  the 
highest  temperature  at  which  a-ferrite  forms  during 
cooling),  but  this  austenite  possesses  low  hardenability. 
On  the  other  hand,  austenitization  becomes  incomplete 
as  one  enters  the  intercritical  region,  within  which  the 
maximum  exposure  temperature  is  between  Acx  and  Ac3. 
However,  the  attendant  austenite  possesses  higher  harde¬ 
nability.  Consequently,  the  volume  fraction  of  the 
freshly  formed  martensite  first  increases  (due  to  the 
dominating  effect  of  the  increased  austenite  hardenabili¬ 
ty),  then  passes  through  a  peak  value  and  finally  de¬ 
creases  (due  to  a  lower  volume  fraction  of  austenite 
available  for  transformation  to  martensite)  as  the  maxi¬ 
mum  exposure  temperature  decreases  from  Ac3  to  Ac\. 
Below  a  maximum  exposure  temperature  corresponding 


to  Ac i ,  no  austenitization  takes  place  during  heating 
and,  therefore,  austenite  —>  (freshly  formed)  martensite 
phase  transformation  takes  place; 

(c)  Variation  of  the  volume  fraction  of  bainite  throughout  the 
FZ  and  within  the  HAZ  (as  a  function  of  distance  from 
the  FZ/HAZ  interface),  Fig.  8(c),  is  merely  a  reflection  of 
the  competition  between  the  austenite  —>  bainite  phase 
transformation  and  the  austenite  — >  (freshly  formed) 
martensite  phase  transformation; 

(d)  Tempered  martensite  does  not  occur  in  the  FZ,  or  within 
the  fully  austenitized  portion  of  the  HAZ,  Fig.  8(e).  In 
the  remainder  of  the  HAZ,  tempered  martensite  volume 
fraction  increases  with  increasing  distance  from  the  FZ/ 
HAZ  interface.  Tempered  martensite  is  the  only  phase 
found  in  regions  of  the  HAZ  with  a  maximum  exposure 
temperature  of  Ac\\  and 

(e)  Moving  from  the  HAZ/unaffected-material  boundary,  the 
grain  size  first  continuously  increases  from  its  initial  value 
(40  pm),  to  a  value  of  ~80  pm  near  the  HAZ/FZ  inter¬ 
face,  Fig.  8(f).  Within  the  FZ,  the  grain  size  first  experi¬ 
ences  a  slight  drop  (the  “chill-zone”  effect)  and  then 
continues  to  increase  towards  the  final  value  of  ~70  pm. 


2.5  Microstructure/Property  Relationship  Module 

As  part  of  the  present  work,  the  microstructure/property 
relationship  module  has  been  upgraded  in  order  to  impart  to  the 
GMAW  process  model  capabilities  for  the  prediction  of 
material  strength  distribution  within  the  weld.  This  upgrade 
utilizes  the  following  three-step  procedure:  (a)  first,  the  key 
contributions  to  the  material  strength  within  different  portions 
of  the  weld  zones  are  identified,  and  the  governing  equations 
for  the  corresponding  strengthening  mechanisms  are  assembled 
and  parameterized;  (b)  next,  a  superposition  scheme  is  devel¬ 
oped  which  will  assess  the  combined  effect  of  different 
strengthening  mechanisms;  and  (c)  finally,  a  functional  rela¬ 
tionship  is  established  between  the  material  strength  and 
hardness  in  order  to  predict  the  spatial  distribution  of  the 
material  hardness  within  the  weld  region.  This  can  then  be  used 
to  validate  the  multi-physics  GMAW  process  model  by 
comparing  the  computed  material  hardness  distribution  with 
its  experimental  counterpart.  Such  validation  is  carried  out  in 
the  portion  of  this  subsection  entitled  “Typical  Results.” 

2.5.1  Main  Contributors  to  the  Weld-Material  Strength. 
To  identify  the  following  main  strength  contributors  in 
prototypical  armor  steels  like  MIL  A46100,  the  crystal 
structure,  chemical  composition,  internal  substructure,  and 
external  morphology  of  the  attendant  phases/micro-constituents 
are  analyzed.  This  analysis  revealed  the  following  main 
strengthening  mechanisms:  (a)  intrinsic  (i.e.,  Peierls  barrier 
control)  strength  of  the  phases/micro-constituents;  (b)  solid- 
solution  strengthening;  (c)  alloy-carbide/cementite  precipitation 
hardening;  (d)  dislocation  substructure-based  strengthening; 
and  (e)  strengthening  by  grain  boundaries  and  phase  interfaces. 
A  brief  overview  of  key  phenomena/processes  pertaining  to 
these  strengthening  mechanisms  and  the  governing  functional 
relationships  will  be  provided  in  the  remainder  of  this  section. 

Intrinsic  Strength  of  the  Phases /Micro-constituents.  Atomic- 
scale  friction  experienced  by  the  gliding  dislocations  during 
plastic  deformation  is  the  main  contributor  to  the  intrinsic 
strength  of  each  phase/micro-constituent  (Ref  15). 
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Fig.  7  Spatial  distribution  of  (a)  equivalent  plastic  strain;  and  (b)  residual  von  Mises  equivalent  stress  over  a  transverse  section  of  the  weld  and 
the  workpiece  region  adjacent  to  the  weld 


Solid-Solution  Strengthening.  The  elastic  (size  and  modulus 
mismatch  type  of)  interactions  between  the  solute  atoms  and 
dislocations  are  assumed  to  be  the  main  controlling  factor  for 
solid-solution  strengthening  (Ref  16).  The  contribution  of  this 
mechanism  to  the  material  strength  is  assessed  using  the 
statistical  approach  of  Labusch  (Ref  17),  which  accounts  for  the 
distribution  of  the  atom/dislocation  interaction  forces  and 
distances.  This  approach  enables  evaluation  of  the  increase  in 
the  critical  resolved  shear  stress,  due  to  the  dislocation/solute- 
atom  interactions,  required  for  dislocation  glide. 

Alloy-Carbide/Cementite  Precipitation  Hardening.  Gliding 
dislocations  experience  increased  resistance  as  a  result  of  their 
interactions  with  precipitates,  which  leads  to  precipitation 
hardening.  The  key  factors  that  control  the  extent  of  precipitate 
hardening  include  precipitate  size,  volume  fraction,  number 
density,  and  chemical/mechanical  properties,  as  well  as  the  state 
of  coherency  at  the  precipitate/matrix  interfaces.  Two  mecha¬ 
nisms,  by  which  gliding  dislocations  generally  overcome 
precipitates,  are  (a)  precipitate  shearing/cutting,  a  process 
which  increases  precipitate/matrix  interfacial  area  and  leads  to 
the  formation  of  interfacial  ledges,  stacking  faults,  and  anti¬ 
phase  boundaries  (in  the  case  of  ordered  crystalline  precipi¬ 
tates);  and  (b)  by  looping  around  the  precipitates,  a  process 
which  leaves  behind  precipitates  surrounded  by  dislocation 
loops  and  increases  the  effective  size  of  the  precipitates. 
Increases  in  the  size  of  the  precipitates  leads  to  increases  in  the 
critical  resolved  shear  stress  required  for  precipitate  shearing. 
However,  for  a  given  volume  fraction  of  precipitate,  the  critical 
resolved  shear  stress  required  for  precipitate-bypass  via  the 
looping  mechanism  decreases  with  an  increase  in  the  precipitate 
size.  Consequently:  (a)  the  extent  of  precipitation  hardening  for 


small  precipitates  is  controlled  by  precipitate  shearing;  (b) 
precipitate  hardening  at  large  precipitate  sizes  is  controlled  by 
dislocation  looping  around  the  precipitates;  and  (c)  there  exists 
a  critical  precipitate  size,  at  which,  the  extent  of  precipitate 
hardening  is  the  highest.  The  degree  of  precipitate  hardening 
associated  with  the  shearing  process  is  dependent  on  the  nature 
and  the  extent  of  phenomena/processes  (e.g.,  precipitate/matrix 
interfacial  energy,  character/magnitude  of  the  stress-free  eigen- 
strain  of  coherent  precipitates,  the  extent  of  precipitate/matrix 
elastic-modulus  mismatch,  the  stacking  fault  energy,  the  anti¬ 
phase  boundary  energy,  etc.)  responsible  for  the  resulting 
increase  in  the  critical  resolved  shear  stress  necessary  for 
dislocation  glide.  The  shearing  mode  of  the  precipitate 
hardening  process  is  expected  to  be  controlled  by  the  elastic- 
modulus  mismatch  between  the  matrix  and  the  precipitates. 
This  is  due  to  the  fact  that  incoherent  alloy-carbide/cementite 
precipitates  are  primarily  found  in  MIL  A46100.  Consequently, 
the  magnitude  of  the  matrix/precipitate  elastic-modulus  mis¬ 
match  is  assumed  to  be  the  key  contributor  to  the  associated 
increase  in  the  critical  resolved  shear  stress  for  dislocation 
glide.  The  critical  resolved  shear  stress  required  for  precipitate- 
bypass  by  the  looping  mechanism  is  quantified  using  the 
Orowan- Ashby  approach  (Ref  18).  This  approach  is  based  on 
the  assessment  of  the  stress  required  to  bend  a  dislocation 
segment  into  a  semi-circular  shape  having  a  radius  equal  to  the 
precipitate  half-distance  along  the  dislocation  line. 

Dislocation  Substructure-Based  Strengthening.  Plastic  relaxa¬ 
tion  takes  place,  both  within  the  growing  martensitic  phase  and 
within  the  neighboring  austenite  phase,  due  to  a  large  shape 
change  (i.e.,  lattice-invariant  shear  strain)  that  accompanies 
austenite  — >  martensite  phase  transformation.  This  results  in  a 
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Fig.  8  Spatial  distribution  of  phase  volume  fractions  for  (a)  allotriomorphic  ferrite,  (b)  Widmanstatten  ferrite,  (c)  bainite,  (d)  freshly  formed 
martensite,  (e)  tempered  martensite,  and  (f)  grain  size  (in  pm)  over  a  transverse  section  of  the  weld  and  the  workpiece  region  adjacent  to  the 
weld 


large  occurrence  of  dislocations  in  freshly  formed  martensite, 
which  act  as  glide-obstacles  to  other  dislocations.  Therefore, 
dislocation  glide  requires  an  increase  in  the  critical  resolved 


shear  stress  and,  hence,  an  increase  in  material  strength.  The 
extent  of  this  strengthening  is  a  function  of  the  transformation 
temperature  since  this  temperature  affects  the  strength  of,  and 
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the  extent  of  plastic  relaxation  in,  the  transforming  austenite.  To 
quantify  this  strengthening  mechanism:  (a)  first,  the  approach 
of  Takahashi  and  Bhadeshia  (Ref  1 9)  is  employed  in  order  to 
establish  a  functional  relationship  between  the  transformation- 
induced  dislocation  density,  p,  and  the  transformation  temper¬ 
ature,  T\  (b)  then,  the  approach  of  Keh  and  Weissman  (Ref  20) 
is  utilized  to  correlate  the  increase  in  critical  resolved  shear 
stress  with  the  increase  in  transformation-induced  dislocation 
density;  and  (c)  finally,  the  approach  of  Leslie  et  al.  (Ref  21)  is 
used  to  assess  the  effect  of  tempering,  which  causes  disloca- 
tion-rearrangement/recovery/partial-annihilation. 

Strengthening  by  Grain  Boundaries  and  Phase  Inter¬ 
faces.  Grain  boundaries  and  phase  interfaces  are  fundamentally 
different  crystalline  defects  in  comparison  to  the  previously 
analyzed  obstacles  to  dislocation  glide  (i.e.,  substitutional/ 
interstitial  alloying  atoms,  precipitates,  and  “forest”  disloca¬ 
tions).  This  difference  stems  from  the  fact  that  grain  boundaries 
and  phase  interfaces  act  as  impenetrable  barriers  to  dislocation 
glide.  Consequently,  during  plastic  deformation,  dislocations 
tend  to  pile-up  at  these  boundaries  and  interfaces.  The 
subsequent  spread  of  plastic  deformation  depends  on  the 
nucleation  and  glide  of  dislocations  in  the  adjacent  grains/ 
phase-particulates  (promoted  by  the  stress-concentration  effects 
arising  from  the  dislocation  pile-ups).  Simple  mechanics 
analysis  of  the  dislocation/boundary  interactions  reveals  that 
the  larger  the  grain/phase-particulate  size,  the  larger  is  the 
dislocation  pile-up  that  can  be  accommodated,  the  larger  are  the 
pile-up-induced  stress-concentration  effects  and,  hence,  the 
lower  is  the  applied  stress  at  which  new  dislocations  are 
nucleated.  In  other  words,  material  strength  decreases  with  an 
increase  in  the  grain/phase-particulate  size.  In  the  microstruc¬ 
ture/property  relationship  module,  this  effect  is  included  using 
the  Hall-Petch  relation  (Ref  22).  Due  to  the  complex  multi¬ 
phase  nature  of  MIL  A46100,  one  can  expect  different  grain 
boundaries  and  phase  interfaces  in  MIL  A46100  (e.g.,  prior- 
austenite  grain  boundaries,  martensite  lath  boundaries,  bainitic 
sheave  boundaries,  and  the  associated  two-phase  interfaces). 
For  all  of  these,  it  is  assumed  that  the  Hall-Petch  relationship  is 
valid. 

2.5.2  Strength  Superposition.  In  this  section,  an  ap¬ 
proach  is  proposed  for  combining  contributions  of  the  strength¬ 
ening  mechanisms,  identified  above,  to  the  total  local  material 
strength  within  the  weld.  The  approach  involves  the  following 
steps: 

(a)  first,  it  is  recognized  that  the  intrinsic  (i.e.,  Peierls  bar¬ 
rier  control)  strengthening  mechanism  provides  the  base¬ 
line  level  of  the  material  strength  while  the  remaining 
four  strengthening  mechanisms  provide  the  respective 
increments  in  the  material  strength  (as  quantified  by 
their  respective  critical  resolved  shear  strengths  for  dislo¬ 
cation  glide); 

(b)  to  obtain  the  total  local  material  strength  within  a  crys¬ 
talline  phase,  the  contributions  of  the  latter  four 
strengthening  mechanisms  are  superposed  and  then 
added  to  the  base-material  strength;  and 

(c)  the  superposition  principle  is  employed  in  the  following 
way:  (i)  for  mechanisms  operating  at  the  same  length 
scale  (e.g.,  solid-solution  strengthening  effects  associated 
with  different  alloying  elements),  the  square  roots  of  the 
corresponding  critical  resolved  shear-strength  increments 


are  first  summed  and  then  squared  (Ref  15);  and  (ii)  for 
mechanisms  operating  at  different  length  scales  (e.g.,  so¬ 
lid-solution  strengthening  and  precipitation  hardening),  a 
simple  linear-superposition  of  the  corresponding  critical 
resolved  shear- strength  increments  is  used. 

The  procedure  described  above  yields  the  total  critical 
resolved  shear  stress  for  dislocation  glide  (i.e.,  the  material 
shear  strength).  To  obtain  the  material  normal  strength,  the 
material  shear  strength  is  multiplied  by  a  factor,  y/3,  in 
accordance  with  the  von  Mises  yield  function. 

It  should  be  noted  that  the  superposition  procedure  described 
above  can  yield  the  strength  of  each  of  the  possible  crystalline 
phases  present  at  a  given  location  within  a  steel  weld.  However, 
the  local  weld  microstructure  typically  does  not  consist  of  a  single 
phase  but  rather  involves  a  number  of  coexisting  phases.  The 
volume  fraction  of  each  of  these  phases  is  predicted  by  the 
microstructure-evolution  module.  To  obtain  the  total  local 
strength  at  a  multi-phase  location  within  the  weld,  a  simple 
rule-of-mixtures  is  employed  within  the  present  module. 

2.5.3  Strength  Versus  Hardness  Relationship.  The 
superposition  methods  described  above  can  be  used  to  predict 
the  normal  strength  at  a  given  material  point  within  the  weld. 
The  resulting  spatial  distribution  of  material  strength  within  the 
weld  can  be  used  to  assess  the  overall  structural  performance  of 
the  weldment.  However,  the  knowledge  of  the  material  strength 
distribution  cannot  be  generally  used  to  validate  the  multi¬ 
physics  GMAW  process  model.  This  validation  would  simply 
require  fabrication  and  testing  of  a  large  number  of  micron¬ 
sized  tensile/compression  specimens.  More  frequently,  it  is  a 
common  practice  to  use  a  distribution  of  indentation  hardness 
(or,  more  precisely,  micro/nano-hardness)  over  the  transverse 
(or  some  other)  section  of  the  weldment,  to  characterize  the 
distribution  of  mechanical  properties  throughout  the  weld 
region.  Since  the  microstructure/property  prediction  module 
yields  the  spatial  distribution  of  the  material  strength  while 
experimental  data  typically  pertain  to  the  spatial  distribution  of 
the  material  hardness,  it  is  necessary  to  establish  the  strength- 
to-hardness  conversion  relation  before  this  portion  of  the 
GMAW  process  model  can  be  validated. 

Ashby  and  Jones  (Ref  23)  suggested  that  material  hardness 
(expressed  in  the  same  units  as  the  normal  strength)  can  be 
calculated  by  multiplying  the  normal  strength  by  a  factor  of  3.0. 
However,  Vickers  micro-hardness  tests  typically  quantify 
material  hardness  in  terms  of  the  so-called  Vickers  hardness 
number  (VHN)  which  is  the  material  hardness  expressed  in 
units  of  kgf/mm2.  To  obtain  VHN,  the  material  hardness, 
expressed  in  MPa,  should  be  divided  by  a  factor  of  9.81. 

2.5.4  Typical  Results.  Figure  9(a)  shows  the  spatial  dis¬ 
tribution  of  the  MIL  A46100  VHN  within  the  weld  region, 
obtained  under  the  same  GMAW  process  and  base-/filler-metal 
conditions  as  the  results  presented  in  Fig.  8(a)-(f).  The  results 
displayed  in  Fig.  9(a)  show  that:  (a)  the  hardness  of  the  material 
within  the  FZ  (the  as-cast  microstructure  of  which  is  dominated 
by  allotriomorphic  and  Widmanstatten  ferrite)  is  lowered 
relative  to  the  hardness  of  the  as-received  material, 
VHN  =  750;  (b)  the  hardness  within  the  (previously  men¬ 
tioned)  “intercritical  region”  is  substantially  higher  than  that  of 
the  as-received  material;  and  (c)  in  general,  both  the  low¬ 
hardness  and  high-hardness  regions  have  to  be  critically 
examined  since  they  lead  to  negative  consequences.  That  is, 
the  low-hardness  region  is  expected  to  compromise  the  load- 
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Fig.  9  Typical  spatial  distribution  of  the  material  VHN  hardness  within  the  weld  region  of  MIL  A46100  as:  (a)  computed  in  the  present  work; 
and  (b)  as  reconstructed  using  the  results  reported  in  Ref  3 1 


bearing  capacity  of  the  weldment  while  the  high-hardness 
region,  through  the  accompanying  loss  of  toughness/ductility, 
may  compromise  weldment  reliability/durability. 

In  the  remainder  of  this  section,  an  attempt  is  made  to 
validate  the  micro  structure/property  relationship  module.  To¬ 
ward  that  end,  the  results  displayed  in  Fig.  9(a)  are  compared 
with  their  experimental  counterparts  (reported  in  Ref  24  and 
recreated)  in  Fig.  9(b). 

A  comparison  of  the  results  depicted  in  Fig.  9(a)  and  (b) 
reveals  that: 

(a)  in  the  HAZ  regions  which  are  located  further  away  from 
the  FZ/HAZ  interface,  the  agreement  between  the  two 
sets  of  results  is  reasonable;  and 

(b)  on  the  other  hand,  some  differences  in  the  two  sets  of 
results  are  observed  in  the  remainder  of  the  weld  region. 
These  differences  can  be  attributed  to  the  fact  that,  while 
in  the  present  work  the  base  and  the  filler-metal  chemis¬ 
tries  are  identical,  in  Ref  24  the  chemical  composition 
of  the  filler-metal  was  somewhat  different  from  that  of 
the  base-metal.  Hence,  within  the  FZ  and  the  region  of 
the  HAZ  adjacent  to  the  FZ  within  which  solid-state  dif¬ 
fusion  was  active,  the  weld  chemistry  was  different  for 
the  results  displayed  in  Fig.  9(a)  and  (b). 


3.  Upgrade  of  the  GMAW  Process  Model 

The  main  objective  of  the  present  work,  as  stated  earlier, 
involves:  (a)  improving  the  predictive  capability  of  the  model 
relative  to  the  mechanical  properties  controlling  the  ballistic- 
limit/penetration-resistance;  and  (b)  imparting  ballistic-limit 
prediction  capabilities  to  the  model.  In  the  remainder  of  this 
section,  these  two  aspects  of  the  GMAW  process  model 
upgrade  are  discussed  separately.  In  the  present  work,  the 
ballistic-limit,  for  a  specific  ballistic-impact  scenario,  is  quan¬ 
tified  using  the  velocity  parameter  V50.  V50  is  defined  as  the 


projectile  incident-velocity  at  which  the  likelihood  of  target 
penetration  is  50%. 

3.1  Upgrade  of  the  Microstructure/Property  Relationship 
Module 

As  mentioned  earlier,  the  microstructure/property  relation¬ 
ship  module  was  originally  capable  of  predicting  strength/ 
hardness  of  different  phases  and  the  overall  material  strength/ 
hardness  distribution  throughout  the  weld.  For  calibrating  and 
validating  the  GMAW  process  model  against  the  available 
experimental  data,  these  predictive  capabilities  are  highly 
critical.  With  respect  to  predicting  the  ballistic-limit  distribution 
throughout  the  weld,  however,  these  capabilities  are  insuffi¬ 
cient.  These  shortcomings  of  the  original  microstructure/ 
property  relationship  module  are  addressed  in  the  present  work 
by  imparting  to  the  model  the  capabilities  of  predicting  the 
weld-spatial  distributions  of  the  following  material  properties: 
(a)  stiffness,  (b)  volume-based  specific  heat,  (c)  exponent  of 
strain-hardening,  (d)  exponent  of  strain-rate  sensitivity,  and  (e) 
material  shear-strength  temperature-dependence.  It  should  be 
noted  that  for  complete  validation  of  the  upgraded  microstruc¬ 
ture,  one  should  validate  the  model  predictions  using  the 
appropriate  experimental  measurements.  This  is  a  formidable 
task  since  the  material  properties  in  question  can  substantially 
vary  throughout  the  weld  region.  Consequently,  experimental 
validation  of  the  present  model  would  require  testing  of  a  large 
number  of  micron-sized  specimens  extracted  from  different 
portions  of  the  weld,  and  was  considered  to  exceed  the  scope  of 
the  present  work.  It  should  be  noted  that  work  is  underway  to 
provide  experimental  validation  of  the  upgraded  microstruc¬ 
ture/property  relationship  module.  However,  in  this  work  a 
simpler  steel  grade  is  being  used  to  reduce  the  extent  of 
experimental  effort  needed. 

The  physical  and  mathematical  meanings  of  the  material 
properties  (c)-(e),  mentioned  above,  will  be  provided  in  the 
next  subsection.  In  addition,  an  attempt  will  be  made  in  the  next 
subsection  to  demonstrate  that  the  spatial  distributions  of  these 
material  properties  must  be  known  in  order  to  predict  the 
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ballistic  limit  and  its  spatial  distribution  throughout  the  weld.  In 
the  remainder  of  this  subsection,  brief  descriptions  will  be 
provided  of  the  procedures  used  to  assess  the  aforementioned 
five  material  properties. 

Material  Stiffness :  The  material  stiffness  properties  of  steels 
(and  their  constituent  phases)  are  typically  considered  as  being 
isotropic  and  microstructure-insensitive.  Following  this 
assumption,  typical  (constant)  values  are  assigned  to  the 
Young’s  modulus  (210  GPa)  and  Poisson’s  ratio  (0.3)  (Ref 

25)  in  the  present  work. 

Volume-Based  Specific  Heat'.  The  volume-based  specific 
heats  of  steels  (and  their  constituent  phases)  are  also  generally 
assumed  to  be  isotropic  and  microstructure-insensitive.  Conse¬ 
quently,  a  typical  (constant)  value  for  the  volume-based  specific 
heat  (3.7  MJ/m3  K)  is  used  (Ref  25)  in  the  present  work. 

Exponent  of  Strain-Hardening'.  The  exponent  of  strain¬ 
hardening  is  generally  considered  to  be  a  micro  structure- 
sensitive  material  property.  In  addition,  this  property  is  assumed 
to  be  inversely  related  to  the  material  strength/hardness  (Ref 

26) .  Figure  10  shows  the  functional  relationship  between  the 
material  yield  strength  and  its  exponent  of  strain-hardening  for 
MIL  A46100  derived  in  the  present  work. 

Exponent  of  Strain-Rate  Sensitivity'.  While  the  exponent  of 
strain-rate  sensitivity  is  generally  found  to  be  moderately 
sensitive  to  the  material  microstructure,  it  is  assumed  to  be 
constant  within  the  weld  (primarily  due  to  a  lack  of  the  relevant 
high-strain-rate  data  for  different  crystalline  phases  and  micro¬ 
constituents  within  the  weld).  The  procedure  recently  proposed 
in  Ref  13  was  utilized  to  determine  a  mean  exponent  of  strain- 
rate  sensitivity  throughout  the  weld,  and  yielded  a  strain-rate 
sensitivity  exponent  of  0.005  for  MIL  A46100  steel.  Since  this 
value  is  relatively  low  in  comparison  to  the  value  of  exponents 
of  strain-rate  sensitivity  found  in  most  metals,  MIL  A46100  can 
be  considered  as  a  weakly  strain-rate-dependent  material. 

Material  Shear-Strength  Temperature-Dependence :  The 
material  shear-strength  temperature-dependence  is  generally 
found  to  be  a  microstructure-dependent  property.  Unfortu¬ 
nately,  there  is  a  lack  of  relevant  experimental  high-temperature 
data  for  this  property  in  different  crystalline  phases  and  micro¬ 
constituents  within  the  weld.  Consequently,  and  following  the 


Fig.  10  Derived  functional  relationship  between  the  material  yield 
strength  and  its  strain-hardening  exponent  for  MIL  A46100 


procedure  outlined  in  Ref  1-3,  this  property  was  assumed  to 
take  on  different  constant  values  within  the  FZ  and  HAZ.  To 
assess  the  mean  values  of  the  material  shear-strength  temper¬ 
ature-dependence  within  the  two  weld-zones,  a  procedure  based 
on  thermally  activated  glide  of  dislocations  within  the  repre¬ 
sentative  microstructures  of  the  two  weld-zones,  originally 
proposed  in  Ref  15,  was  employed  in  the  present  work.  Within 
the  FZ  and  HAZ  of  MIL  A46100,  this  procedure  yielded  the 
(room  temperature)  shear-strength  temperature-dependence 
values  of  —0.17  and  —0.15  MPa/K,  respectively.  The  lower 
magnitude  of  the  shear- strength  temperature-dependence  in  the 
HAZ  can  be  attributed  to  the  operation  of  auto-tempering 
processes  within  this  region.  As  far  as  the  (room  temperature) 
shear  strength  of  the  subject  material  MIL  A46100  is  con¬ 
cerned,  its  as -received  value  of  577  MPa  has  been  taken  from 
Ref  4.  On  the  other  hand,  within  the  weld,  the  results  obtained 
in  the  present  work  revealed  that  the  (room  temperature)  shear 
strength  varies  in  an  approximately  500-650  MPa  range. 

3.2  Creation  of  the  Ballistic-Limit  Prediction  Module 

The  procedure  used  to  construct  the  sixth  (ballistic-limit 
prediction)  module  of  the  GMAW  process  model  includes  the 
following  steps:  (a)  identification  of  the  failure  mechanisms 
controlling  the  ballistic  limit;  (b)  for  each  of  the  failure 
mechanisms  identified,  a  physically  based  model  is  developed 
to  predict  the  associated  ballistic  limit;  (c)  the  associated 
governing  equations  and  their  parameterizations  are  next 
incorporated  into  the  sixth  module;  and  (d)  last,  the  sixth 
module  is  utilized  to  predict  the  spatial  distribution  of  the 
ballistic  limit  within  a  GMAW  joint  in  MIL  A46100. 

3.2.1  Identification  of  the  Ballistic-Limit-Controlling 
Failure  Mechanisms.  Based  on  a  review  of  the  public- 
domain  literature  conducted  as  part  of  the  present  work,  the 
following  main  failure  modes  controlling  the  ballistic  limit  in 
metallic  welds  were  identified  (Ref  27-30):  (a)  front/impacted 
and  back/rear  face  petaling,  Fig.  11(a)  and  (b);  (b)  ductile  hole- 
enlargement,  Fig.  12;  (c)  plugging,  Fig.  13;  and  (d)  spalling, 
Fig.  14(a)  and  (b).  Both  the  weld  factors  (e.g.,  weld  soundness/ 
quality,  relative  and  absolute  sizes  of  the  FZ  and  HAZ,  weld 
micro  structure,  spatial  distribution  and  type  of  micro  structure 
within  the  weld,  etc.)  and  the  impactor  properties  (e.g., 
projectile’s  mass,  shape,  size  (relative  to  the  target  thickness), 
and  hardness  (relative  to  that  of  the  target-material))  have  been 
found  to  control  the  dominant  failure  mode  in  a  specific 
ballistic-impact  scenario.  Since  detailed  descriptions  of  the  key 
features  of  the  aforementioned  four  failure  mechanisms  can  be 
found  in  Ref  27-30,  the  same  details  will  not  be  presented  here. 

3.2.2  Modeling  of  the  Failure  Mechanisms.  This  sub¬ 
section  focuses  on  assembling  the  relations  quantifying  the 
contribution  of  the  four  failure  modes  to  the  target  ballistic 
limit. 

Front/Impacted  and  Back/Rear  Face  Petaling'.  Since  this 
mode  of  failure  makes  a  small  overall  contribution  to  the 
ballistic  limit  of  thick  armor-plates  (Ref  31),  the  role  of  this 
failure  mechanism  in  the  overall  ballistic  performance  of 
GMAW-joints  was  not  modeled  and/or  assessed  in  the  present 
work. 

Ductile  Hole-Enlargement'.  To  model  this  failure  mecha¬ 
nism,  the  so-called  “Cavity  Expansion  Theory  (CET)”  (Ref  29, 
30,  32-34)  is  employed  here.  This  theory  treats  target  penetra¬ 
tion  by  a  (rigid,  ogive-nosed)  projectile  as  a  problem  of  radial 
expansion  of  the  associated  (cylindrically  shaped)  cavity  within 
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the  target.  To  derive  the  ballistic-limit  relation,  the  following 
three-step  procedure  is  employed:  (a)  the  governing  mass  and 
momentum  conservation  equations  are  first  combined  with  the 
elastic-plastic  material  constitutive  relations  for  the  target,  to 
establish  a  functional  relation  between  the  axial  force  opposing 
target  penetration  and  the  cavity-expanding  radial  stress  (as); 
(b)  the  latter  quantity,  as,  is  next  correlated  with  the  target- 
material  stiffness  and  strength  (including  strain-hardening) 
parameters  and  with  the  cavity  expansion  rate;  and  (c)  by 
employing  Newton’s  second  law,  and  by  ignoring  the  contri¬ 
bution  of  the  radial  expansion  inertia  to  <3S ,  the  ballistic  limit 
V50  is  defined  as  the  projectile  minimal  incident-velocity 
required  for  the  projectile  to  traverse  the  entire  target  thickness, 
and  expressed  as  (Ref  35): 


^27 za2has^ 1 ^ 


(Eql) 


where  a  is  the  projectile  shank  radius,  m  the  projectile  mass, 
and  h  the  target  thickness. 

For  a  nonlinear  strain-hardening  target-material  with  the 
post-yield  uniaxial  stress,  gu  versus  uniaxial  strain  su  relation 
being  defined  as: 


(Eq  2) 


Fig.  11  Schematic  representation  of  (a)  front  face;  and  (b)  back 
face  petaling  failure  mechanism 


where  Y  is  the  quasi-static  yield  strength,  E  is  Young’s  modu¬ 
lus,  and  n  is  the  strain-hardening  exponent,  gs  can  be  ex¬ 
pressed  as 


Fig.  12  Schematic  representation  of  the  ductile  hole-enlargement 
failure  mechanism 


Fig.  13  Schematic  representation  of  the  plugging-failure  mecha¬ 
nism 


(-In  x)n 
(1  -x) 


b  = 


V3Y 

E 

(Eq  3) 


Equations  1  and  3  are  then  combined  to  obtain: 


Equation  4  reveals,  in  accordance  with  experimental  find¬ 
ings,  a  weak  nonlinear  relationship  between  material  strength 
(hardness)  and  the  ballistic  limit  (E5o).  While  one  should,  in 
principle,  use  high- strain-rate  material  properties  to  compute 
V50  in  Eq  4,  it  is  generally  found  that  the  use  of  quasi-static 
material  properties  still  yields  reasonable  predictions  for  V50 
through  Eq  4.  This  observation  is  generally  rationalized  using 
the  fact  that,  under  ballistic  loading,  material  deformation  is 
nearly  adiabatic  (i.e.,  heat  conduction  is  practically  absent),  and 
hence,  ductile  hole-enlargement  is  controlled  by  the  plastic 
deformation  of  the  surrounding  material  at  (practically)  ambient 
temperature  (at  which  high-strength  steels  typically  show  very 
little  strain-rate  sensitivity). 

Plugging :  Shear  cracking  is  generally  assumed  to  cause  this 
mode  of  failure.  Shear  cracking  is  believed  to  be  induced  by  the 
formation  of  adiabatic  plastic-shear  bands  within  the  target- 
plate  material  surrounding  the  advancing  blunt  projectile 
(Ref  27,  29,  36).  To  model  plugging-failure,  the  conservation 
equations  and  material  constitutive  relations  are  combined  with 
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Fig.  14  Schematic  representation  of  the  spalling  failure  mechanism  within  the  (a)  HAZ;  and  (b)  FZ  regions 


one  of  the  adiabatic  plastic-shear  instability  theories  (e.g., 
Ref  36)  in  the  present  work. 

Over  the  last  30  years,  the  phenomenon  of  adiabatic  plastic- 
shear  instability  and  its  subsequent  localization  has  been 
investigated  extensively  Nowadays,  it  is  believed  that  the 
essential  physics  behind  adiabatic  plastic-shear  instability/ 
localization  is  well  understood.  However,  reliable  quantitative 
relations  between  the  onset  of  this  failure  mechanism  and  the 
material  properties  and  loading  conditions  are  required.  The 
onset  of  plastic-shear  strain  instability  is  generally  associated 
with  a  condition  when,  in  the  course  of  loading,  strain-induced 
hardening  becomes  overmatched  by  the  strain-induced  soften¬ 
ing  processes.  Under  such  conditions,  the  local  rate  of  change 
of  the  maximum  shear  stress,  x,  with  the  associated  shear  strain, 
y,  becomes  non-positive,  i.e., 


n 

(~\  dx\ 

m  dy\ 

C  dT  ly,y 

j  dy) 

(Eq  6) 


where  yc  is  a  (minimum)  critical  value  of  the  shear  strain  at 
the  onset  of  plastic-shear  strain  instability  at  which  ^  =  0.  In 
the  high- strain-rate  regime  (  y>100s-1)  and  for  a  typical 
range  of  values  of  the  strain-rate  sensitivity  exponent 
(0  <  m  <  1.0),  Eq  6  can  be  further  simplified  as: 


(Eq  7) 


Equation  7  reveals  that  the  onset  of  plastic-shear  strain 
instability  and  localization  is  delayed  in  materials  in  which: 
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(Eq  5) 


where  a  raised  dot  is  used  to  denote  a  time  derivative  and  T 
denotes  temperature. 

Under  adiabatic  conditions,  the  term  ^  in  Eq  5  can  be  set 
equal  to  x/C,  where  C  is  the  volumetric  specific  heat.  As  far  as 
the  remaining  derivatives  on  the  right-hand  side  of  Eq  5  are 
concerned,  their  functional  relationship  depends  on  the 
x(y,y,  T)  relation  used.  For  instance,  if  strain-hardening  and 
strain-rate  sensitivity  relations  are  both  assumed  to  be  of  a 
parabolic  character  with  the  corresponding  exponents  denoted, 
respectively,  as  n  and  m,  Eq  5  can  be  re-written  as: 


(a)  volumetric  specific  heat  is  relatively  large  and,  hence, 
the  associated  adiabatic  temperature  increase  is  relatively 
small; 

(b)  strain-hardening  effects  are  strong  and  can  effectively 
counteract  the  effect  of  various  softening  processes;  and 

(c)  temperature  sensitivity  of  the  material  strength  is  not 
pronounced,  so  that  local  heating  does  not  significantly 
reduce  the  material  resistance  toward  plastic  deforma¬ 
tion. 

Examination  of  Eq  6  and  7  reveals  that  they  provide 
functional  relationships  between  the  critical  shear  strain  at  the 
onset  of  plastic-shear-strain  localization  and  the  material 
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constitutive  response  (properties).  However,  these  equations  do 
not  explicitly  address  the  target-material  ballistic  limit. 

While  it  is  well-established  that  plugging-failure  is  domi¬ 
nated  by  shear-strain  localization,  the  task  of  establishing 
functional  relationships  between  the  associated  ballistic  limit 
and  the  material’s  propensity  to  plastic- shear-instability  (as 
quantified  by  the  critical  shear  strain,  yc)  is  formidable.  For 
example,  while  an  increase  in  material  hardness  is  generally 
found  to  increase  the  plugging-failure-controlled  ballistic  limit, 
this  functional  relationship  is  not  reflected  in  Eq  6  and  7. 
Consequently,  one  would  expect  that  high-hardness  and  high- 
strain-hardening  target-materials  would,  in  general,  possess 
superior  ballistic  performance.  However,  in  the  case  of  welded 
joints,  there  is  a  trade-off  between  material  hardness  and  its 
strain-hardening  potential  (i.e.,  these  two  material  properties  are 
inversely  related  and,  hence,  cannot  be  changed  indepen¬ 
dently).  Furthermore,  it  should  be  noted  that  the  ballistic-limit 
scales  with  the  energy  absorbed  by  the  target,  which  can  be 
defined  as  an  integral  of  stress  over  the  associated  strain  path 
until  localization.  Consequently,  at  the  lowest  hardness  levels 
and,  thus,  the  highest  strain-hardening  levels,  the  onset  of 
plastic-shear-strain  instability  is  greatly  delayed.  However,  the 
overall  projectile  kinetic  energy  absorption  (i.e.,  the  ballistic 
limit)  remains  low.  As  the  target-material  hardness  increases,  its 
ballistic  limit  increases  due  to  the  accompanying  higher  ability 
for  energy  absorption.  In  other  words,  the  effect  of  an  increase 
in  the  material  strength/hardness  more  than  compensates  for  the 
accompanying  reduction  in  yc  (caused  by  the  associated 
reduction  in  the  strain-hardening  potential).  At  the  highest 
hardness  levels  and,  thus,  the  lowest  strain-hardening  levels,  the 
target-material’s  ability  to  absorb  projectile  energy  decreases 
due  to  the  fact  that  the  material  becomes  highly  prone  to  plastic 
strain  localization  (a  very  small  value  of  yc). 

The  aforementioned  analysis  reveals  a  complicated  interplay 
between  material  hardness/strength,  its  strain-hardening  behav¬ 
ior  and  its  plugging-failure-controlled  ballistic  limit.  This 
appears  to  be  the  main  reason  for  the  lack  of  a  broadly 
accepted  functional  relationship  between  the  plugging-failure- 
controlled  ballistic  limit  and  the  target-material  constitutive 
response/properties.  To  overcome  this  problem,  a  simple 
phenomenological  model  for  the  plugging-failure-controlled 
ballistic  limit  is  developed  in  the  remainder  of  this  section.  The 
model  is  based  on  the  following  simplifying  assumptions/ 
approximations: 

(a)  the  projectile  (of  a  right-circular  cylindrical  shape)  re¬ 
mains  undeformed  during  impact;  and 

(b)  the  kinetic  energy  of  the  projectile  at  the  onset  of  plug¬ 
ging  is  exactly  balanced  by  the  work  of  plastic  deforma¬ 
tion  in  the  shear-localization  region  surrounding  the 
projectile  as: 

Yc 

2  PProj^Proj^Proj  vh  —  2  7li?proj  A  Target  bRshear  ^  'tdy  •>  (Eq  8) 

0 

where  p  denotes  density,  R  the  radius,  h  the  thickness/height, 
5 R  the  radial  thickness  and  subscripts  “Proj”,  “Target,”  and 
“Shear”  represent  the  appropriate  quantity  in  the  projectile, 
target,  and  the  plastic-shear-localization  zone. 

It  should  be  noted  that  the  present  model  contains  only  one 
adjustable  parameter,  5i?Shear-  This  parameter  has  been  deter¬ 
mined  by  curve-fitting  the  plugging-failure-controlled  ballistic- 


limit  experimental  data  for  MIL  A46100.  Specifically,  V50 
versus  /zTarget  data  are  taken  from  Ref  4.  The  data  used  pertain 
to  the  case  of  a  20  mm  caliber,  53  g,  blunt,  chamfered  right- 
circular  cylindrical  steel  fragment  simulating  projectile  or  FSP 
(described,  in  greater  detail,  in  the  later  portion  of  this  section). 
Thus,  all  the  quantities  with  subscript  Proj  appearing  in  Eq  8 
are  known.  Then,  for  each  pair  of  values  of  V50  and  /^Target? 

shear  is  determined  from  Eq  8.  Finally,  an  average  value  of 
d^shear  =  0.4mm  is  obtained. 

Spalling :  This  mode  of  failure  possesses  the  following  main 
characteristics:  (a)  typically  occurs  at  the  target  back  face;  (b)  is 
associated  with  a  loss  of  the  material  through-the-thickness 
toughness;  (c)  usually  gives  rise  to  a  substantial  increase  in  the 
size  of  the  projectile-exit  hole  (in  the  case  of  target  defeat);  and 
(d)  gives  rise  to  a  loss  in  the  target-material  ballistic  limit 
(particularly  in  the  weld  region  characterized  by  large  hardness 
and,  thus,  low  toughness  levels). 

Although,  over  the  last  30  years,  spall  failure  has  been 
investigated  extensively,  many  aspects  of  this  important  phe¬ 
nomenon  remain  unresolved  (e.g.,  Ref  37).  This  is  the  main 
reason  that  presently,  there  is  no  broadly  accepted  functional 
relation  for  the  spall-failure-dominated  ballistic  limit.  In  the 
present  work,  no  attempt  will  be  made  to  derive  this  functional 
relationship.  This  decision  was  based  on  the  following  reasons: 

(a)  The  overall  contribution  of  the  spalling  failure  to  the 
ballistic  limit  of  metallic  welds  is  relatively  small  in 
comparison  to  the  contributions  associated  with  the  duc¬ 
tile  hole-enlargement  and  plugging-failure  mechanisms; 
and 

(b)  Modeling  of  the  spalling-controlled  ballistic  limit  re¬ 
quires  knowledge  of  the  spatial  distribution  of  material 
fracture  toughness  throughout  the  weld.  However,  the 
micro  structure/property  relationship  module,  even  after 
the  inclusions  of  the  upgrades  developed  in  the  present 
work,  is  unable  to  provide  predictions  regarding  the  spa¬ 
tial  distribution  of  fracture  toughness  within  the  weld. 

3.2.3  Typical  Results.  The  discussion  in  the  previous 
subsection  clearly  established  that  the  geometry  of  the  projec¬ 
tile  affects  whether  the  ballistic  limit  of  the  target  is  controlled 
by  hole-enlargement  failure  or  plugging-failure.  For  this  reason, 
separate  predictions  regarding  the  spatial  distribution  of  v50 
within  the  MIL  A46100  GMAW  joint  are  made  for  these  two 
modes  of  target  penetration. 

Ductile  hole-enlargement  failure  mode :  In  this  case,  the  V50 
velocity  is  calculated  using  Eq  4  (along  with  the  local  values 
for  the  material  Young’s  modulus,  yield  strength,  and  exponent 
of  strain-hardening).  Only  one  type  of  projectile,  a  7.62  mm 
caliber,  8.2  g,  ogival-nose,  steel-jacketed,  hard  tungsten  core 
(1400  HV),  and  armor  piercing  (AP)  projectile,  was  used  in  the 
analysis.  The  Young’s  modulus  of  the  target  is  treated  as  a 
microstructure-insensitive  property  and,  hence,  assigned  a 
constant  value  throughout  the  weld.  The  results  displayed  in 
Fig.  9(a)  are  used  to  assign  the  spatial  distribution  of  the  target- 
material  yield  strength/hardness  within  the  joint.  The  local 
value  of  the  exponents  of  strain-hardening  throughout  the  weld 
was  obtained  by  combining  the  results  displayed  in  Fig.  9(a) 
with  the  functional  relationship  depicted  in  Fig.  10. 

The  spatial  distribution  of  V50  yielded  by  this  procedure  is 
displayed  in  Fig.  15(a).  The  results  displayed  in  this  figure  were 
obtained  for  the  case  of  a  12.7  mm  thick  MIL  A46100 
armor-plate,  welded  under  the  same  GMAW  process  and 
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(b) 


Fig.  15  Typical  spatial  distributions  of  the  ^50  (m/s)  for  the  case  of  (a)  ductile  hole-enlargement  failure;  and  (b)  plugging-failure-controlled  tar¬ 
get  penetration/ballistic-limit 


base-/filler-metal  conditions  as  the  results  presented  in 
Fig.  8(a)-(f).  Examination  of  Fig.  15(a)  reveals  that: 

(a)  there  is  a  region  within  the  FZ  within  which  the  ballistic 
limit  has  been  reduced  (by  ca.  15%)  below  the  level  ob¬ 
served  in  the  base-metal  ca.  1455  m/s; 

(b)  there  is  a  “critical  region”  within  the  HAZ  within  which 
V50  has  dropped  by  7-8%  relative  to  its  base-metal  level. 
The  presence  of  such  critical  regions  within  the  HAZ 
has  been  reported  in  Ref  24;  and 

(c)  in  the  particular  case  investigated  here,  the  region  asso¬ 
ciated  with  a  lowered  value  of  V50  represents  a  major 
portion  (>10%)  of  the  weld.  Therefore,  one  can  expect 
that  the  presence  of  the  GMAW  joint  will,  in  the  case  of 
direct  impact,  compromise  the  ballistic  performance  of 
the  armor-plate. 

Plugging-failure  mode :  In  this  case,  the  V50  velocity  is 
calculated  using  Eq  7  and  8  (along  with  the  local  values  of  the 
volume-based  heat  capacity,  material  yield  strength,  exponent 
of  strain-hardening,  and  the  material  shear- strength  tempera¬ 
ture-dependence  ff;|Yp.  A  single  projectile  case,  a  20  mm 
caliber,  53  g,  blunt,  and  chamfered  right-circular  cylindrical 
steel  FSP,  is  used  in  the  analysis.  As  mentioned  earlier,  the 
volume-based  heat  capacity  of  the  target  is  treated  as  a 
microstructure-insensitive  property  and  assigned  a  constant 
value  throughout  the  weld.  The  spatial  distribution  of  the  target- 
material  yield  strength  within  the  joint  is  again  taken  from 
Fig.  9(a)  and  the  exponent  of  strain-hardening  is  obtained  by 
combining  the  results  displayed  in  Fig.  9(a)  with  the  functional 
relationship  depicted  in  Fig.  10.  As  far  as  the  temperature 
sensitivity  of  the  material  strength  is  concerned,  it  is  assumed, 
as  mentioned  earlier,  to  have  constant  but  different  values 
within  the  FZ  and  HAZ  (Ref  38-41). 

The  procedure  described  above  yielded  the  spatial  distribu¬ 
tion  of  V5Q  as  shown  in  Fig.  15(b).  These  results  were  obtained 
for  the  case  of  a  12.7  mm  thick  MIL  A46100  armor-plate, 


welded  under  the  same  GMAW  process  and  base-/filler-metal 
conditions  as  the  results  presented  in  Fig.  8(a)-(f).  Examination 
of  the  results  displayed  in  Fig.  15(b),  and  their  comparison  with 
the  results  displayed  in  Fig.  15(a),  reveals  that: 

(a)  the  overall  spatial  distribution  of  the  ballistic  limits  with 
respect  to  the  plastic-shear-localization  and  ductile  hole- 
enlargement  failure  mechanisms  are  fairly  similar;  and 

(b)  the  F50  values  pertaining  to  the  plastic-shear  localization 
failure  mode  are  consistently  lower  than  their  ductile 
hole-enlargement  failure  counterparts.  This  finding  is 
fully  consistent  with  the  earlier  discussion  which  indi¬ 
cated  that,  at  high-strength  levels,  steels  become  more 
susceptible  to  shear-localization-induced  failure  under 
ballistic-impact  conditions. 


4.  Summary  and  Conclusions 

Based  on  the  work  presented  and  discussed  in  the  present 
manuscript,  the  following  main  summary  remarks  and  conclu¬ 
sions  can  be  made: 

(1)  The  five-module  multi-physics  computational  model  for 
the  conventional  GMAW  joining  process,  developed  in 
our  recent  work,  has  been  expanded  to  include  a  sixth 
module.  Within  this  module,  the  model  is  imparted  the 
capabilities  for  predicting  the  spatial  distribution  of  the 
mechanical  properties  controlling  the  ballistic  limit  (i.e., 
penetration-resistance)  of  the  weld. 

(2)  The  new  module  enables  prediction  of  the  all-metal 
welded-armor  ballistic-limit  under  the  condition  that  this 
limit  is  controlled  by  either  ductile  hole-enlargement  or 
plugging-failure  mechanisms. 

(3)  To  demonstrate  the  utility  of  the  upgraded  GMAW 
process  model,  the  case  of  a  prototypical  low-alloy, 
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high-hardness  armor-grade  martensitic  steel,  MIL  A46100, 
welded-armor  (with  the  filler-metal  used  being  also  MIL 
A46100)  is  analyzed. 

(4)  The  work  presented  in  this  manuscript  clearly  revealed 
the  capacity  of  the  multi-physics  GMAW  process  model 
to  establish  functional  relationships  between:  (a)  the  pro¬ 
cess  parameters,  e.g.,  input  welding  voltage,  electrode- 
to-workpiece  distance,  filler-metal  feed-rate,  welding- 
gun  speed,  etc.;  (b)  the  resulting  spatial  distribution  of 
the  material  micro  structure  (as  represented  by  the  vol¬ 
ume  fractions  of  various  crystallographic  phases  and  mi¬ 
cro-constituents  as  well  as  by  the  prior-austenite  grain 
size);  (c)  the  associated  mechanical  properties  (e.g., 
strength/hardness,  stiffness,  exponent  of  strain-harden¬ 
ing,  etc.);  and  (d)  the  ballistic-limit/penetration-resis- 
tance  within  the  fusion  and  the  heat-affected  zones  of 
the  weld. 
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